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Abstract 

We investigate the dark matter prospects of supersymmctric models with nonuniversal gaug- 
ino masses. We find that for very particular values of the ratio of soft supersymmetry-brcaking 
gaugino masses, M2/M1, an enhanced coannihilation efficiency between the lightest chargino 
and the lightest neutralino occurs, allowing for scalars with masses well above the normally 
accepted limits for viable dark matter in the universal case. As a specific example, we investi- 
gate models of hidden sector gaugino condensation. These models exhibit high scalar masses, 
previously thought dangerous, and the requisite freedom in the ratio of gaugino masses. The cos- 
mologically viable regions of parameter space are investigated, allowing very specific statements 
to be made about the content of the supersymmetry-breaking hidden sector. 
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It has long been held that one of the prime virtues of supersymmetry as an explanation of 
the hierarchy problem is that it tends to also provide a solution to the dark matter problem as 
a nearly automatic consequence of R-parity conservation. The lightest supersymmetric particle 
(LSP) is then stable and, as it tends to be a neutral gaugino, it will typically have the right mass 
and annihilation rate in the early universe to provide sufficient mass density today to account for 
observations suggesting p tot ~ p ciit |]] . 

This paper initially investigates the dark matter implications of the most widely studied bench- 
mark in supersymmetric phenomenology, the constrained Minimal Supersymmetric Standard Model 
(cMSSM). We emphasize that the cMSSM fails to solve the dark matter problem over most of its 
parameter space with the exception of certain very special patterns of soft supersymmetry-breaking 
terms. These patterns must constrain the neutralino (a fermion) to have a very specific mass re- 
lationship to an unrelated boson such as the lightest Higgs or the stau. Barring these fine-tuned 
relationships, the cMSSM predicts too much dark matter - thus bringing it into conflict with direct 
measurements of the age of the universe ||, |3| . We find this failure to be due, in part, to the cMSSM 
constraint on the gaugino mass ratio M<ijM\. Eliminating this assumption of universal gaugino 
masses uncovers new regions of parameter space that allow for cosmologically allowed, and often 
experimentally preferred, values of the neutralino relic density. 

This paper is organized as follows: in Section p] we present our methodology in the context of the 
cMSSM with its standard minimal supergravity (mSUGRA) universal soft supersymmetry-breaking 
terms. Subsequently, in Section ^, we relax our assumption of universal gaugino masses. While the 
relic density implications of nonuniversal gaugino masses, and in particular the role of M%IM\ in 
dark matter phenomenology, have been explored previously these past studies have either focused 
on specific models or have not included the important effects of coannihilation between the LSP 
and the lightest chargino [[|, ||].[] 

In Section |3] we consider a specific class of supergravity models derived from heterotic string 
theory which implement supersymmetry breaking through gaugino condensation in a hidden sec- 
tor @, ||, [| as an example of how the general results of Section || can be applied on a model-by- 
model basis. Requiring a cosmologically relevant thermal LSP relic density will imply very specific 
conclusions about the content of the hidden sector of these models. Finally, we conclude and remark 
upon possible extensions of this work. 

1 A noteworthy exception is Ref [H] though it focuses primarily on a purely wino-like LSP scenario with LSP masses 
below Mw 
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1 Universal Gaugino Masses 



The phenomenological consequences of the cMSSM have been studied extensively [jD] ], including the 



cosmological implications of its (presumed stable) LSP [12|. In such a regime the entire low energy 
phenomenology is specified by five parameters: a common gaugino mass M1/2, a common scalar 
mass Mq, a common trilinear scalar A-term Aq, the value of tan (3 and the sign of the //-parameter 
in the scalar potential. These values are defined at some high energy scale, typically taken to be 
the scale of gauge coupling unification Auv ~ 2 x 10 16 GeV. 

To obtain the superpartner spectrum at the electroweak scale the renormalization group equa- 
tions (RGEs) are run from the boundary scale to the electroweak scale |l3|]. In the following, all 
gauge and Yukawa couplings as well as A-terms were run with one loop RGEs while scalar masses 
and gaugino masses were run at two loops to capture the possible effects of heavy scalars and a 
heavy gluino on the evolution of third generation squarks and sleptons. We chose to keep only the 
top, bottom and tau Yukawas and the corresponding A-terms. 

At the electroweak scale Aew = Mz the one loop corrected effective potential Vi_i 00 p = Vtree + 
AVrad i s computed and the effective //-term p, is calculated 



( m H d + Sm H d ) ~ ( m H u + 6m H u ) tan /5 



tan 2 0-1 



-Ml 



(1) 



In equation ([!]) the quantities 5mu u and 5m n d are the second derivatives of the radiative corrections 
AVrad with respect to the up-type and down-type Higgs scalar fields, respectively. These corrections 
include the effects of all third generation particles. If the right hand side of equation is positive 
then there exists some initial value of \i at the high energy scale which results in correct electroweak 
symmetry breaking with Mz = 91.187 GeV. 

The neutralino states and their masses are calculated using the neutralino mass matrix 
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where Mi is the mass of the hypercharge U(l) gaugino at the electroweak scale and M2 is the mass 
of the SU(2) gauginos at the electroweak scale. The matrix (0) is given in the (B, W, H^, H®) basis, 
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where B represents the bino, W represents the neutral wino and H® and H® are the down-type 
and up- type Higgsinos, respectively.^ 

The lightest eigenvalue of this matrix is then typically the LSP and it is overwhelmingly bino- 
like in content over most of the parameter space when tan (3 is low. This is because the cMSSM 
universality constraint on gaugino masses at the high scale of the theory implies M\ ~ \M 2 when 
the masses are evolved to the electroweak scale via the RGEs. Provided |Mi|, \M 2 \ <C which 
is the case for low tan/3, the LSP mass is then dominated by M\ and has a typical bino content 
of <; 99%. We will restrict ourselves to this low tan/3 regime and adopt a value of tan/3 = 3 for 
the remainder of the paper. The dark matter prospects of the high tan (3 limit have been studied 
recently by Feng et al. [fl5|| . 

More generally the content of the LSP can be parametrized by writing the lightest neutralino 

as: 

X 1 =N 11 B + N 12 W + N 13 H% + NuH°, (3) 

which is normalized to + iVj^ + ^13 + ^Yu = 1- Thus by saying that the bino content of the 
lightest neutralino is high, we mean N\\ — 1. 

Given the particle spectrum we compute the thermal relic LSP density with a modified version 



of the software package neutdriver [16]. This program includes all of the annihilation processes 
computed by Drees and Nojiri [jl?]] which are used to compute a thermally averaged cross section 
< (TV >ann arid freeze-out temperature xf = Tp/m x o. From knowledge of the mass of the lightest 
neutralino m x o (assumed to be the LSP), < av > an n and Tp, a relic abundance can be computed 
using the standard approximation 



2 1.07 x 10 9 xp 

1/2 

gj Mpi (a e ff + 3 (b eS - a eS /A) /x F ) 



^ 2 = 1/2 „ ; GeV- 1 , (4) 



where we have expressed the thermally averaged annihilation cross section as an expansion in 
powers of the relative velocity: 

< crv > ann = a cS + b cfi v 2 + ... (5) 

A proper determination of relic LSP densities requires that the above computation be amended 
to include the possible effects of coannihilation fll8||. This occurs when another particle is only 



2 Loop corrections at next-to-leading order jTl|] to this mass matrix have been incorporated and found to have 
little effect on the results that follow. 
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slightly heavier than the lightest neutralino so that both particles freeze out of equilibrium at 
approximately the same temperature. The neutralino can now not only deplete its relic abun- 
dance through annihilation processes such as XiX? — * e + e~ , but also through interactions with the 
coannihilator such as X1X1 ~^ e±z/ e- The extreme importance of including relevant coannihilation 
channels has recently been emphasized for the case of the cMSSM Jl9| , [20| , and in that spirit we 
have added a number of coannihilation channels which are relevant for both the universal gaugino 
mass case [20| as well as the case of nonuniversal gaugino masses to be considered in Section |2[ 

The program neutdriver includes x^Xi coannihilation to W^j [Q] and two (massless) fermions 
//' which we recalculated to account for non-zero fermion masses. We have also included a calcu- 
lation of the process x^Xi ~ * W^Z, and found this channel to often dominate when kinematically 
accessible. Additionally, we have inserted the results of [19] for Xi^ coannihilation to Zt, 77", hr, 
and Hr final states. 

The region of cMSSM parameter space that gives rise to acceptable levels of bino-like LSP relic 
density is given in Figure p] where we have plotted contours of ^ x h 2 = 0.1, 0.3 and ^ x h 2 = 1. Here 
fi x is the fractional LSP matter density relative to the critical density and h is the reduced Hubble 



parameter: h ~ 0.65 [21]. We have chosen M1/2 and Mo as free parameters in the manner of Ellis, 
Falk & Olive [^] with tan (5 = 3, Aq = and positive //-term.0 

This plot is similar to the ones presented in [ 20 1 and we reproduce it here to draw attention to two 
important facts. First, observations suggest that the preferred values for cold dark matter densities 
are in the range 0.1 < f^h 2 < 0.3 |22(] , though including the latest evidence of a cosmological 
constant |J may shift this to 0.06 ^ ^x^ 2 ^ 0-2- This experimental data points to a region of 
the cMSSM parameter space in which both the universal scalar mass and the universal gaugino 
mass are small, on the order of 200 GeV for each (see, for example, the region around point A 
in Figure |l]). In fact, heavy scalars can only be accommodated cosmologically if nature was kind 
enough to arrange the soft-supersymmetry breaking parameters so that the mass of the LSP is 
about half of the mass of the lightest Higgs boson. In that case, though t-channel scalar fermion 
exchange may be suppressed and the annihilation rate into two fermions is not sufficient to deplete 
the LSP density adequately, the annihilation through resonant s-channel exchange of the lightest 
Higgs is efficient enough to provide an acceptable dark matter region insensitive to the scalar mass. 

We have allowed our Mq value to range to as much as a TeV, in contrast to pCf], to accentuate 



3 In our conventions this is the sign of fj, least constrained by the measurement of the branching ratio for b — > 5-7 
events. This is the opposite convention used by the neutdriver package. 
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the point that most of the cMSSM parameter space considered "natural" in the literature predicts 
far too much neutralino relic density. Every point in Figure |] above the ^ x h 2 = 0-3 li ne is already 
excluded by astrophysical measurements of the dark matter density. This rules out almost all of the 
plotted parameter space except for the low-mass region including point A and the thin strip right 
above the region excluded due to a charged LSP (the region including point B). Once we begin to 
impose constraints arising from Higgs searches at LEP, the preferred low mass region in Figure ffi 
begins to be ruled out. The only region which then remains in the {Mq, Mi/ 2 } plane for tan/3 = 3 
is that which is due to coannihilation between the LSP and the lightest stau. Being in this region 
requires a very specific relationship between the gaugino mass parameter and the unrelated scalar 
mass parameter. 

To better illustrate the physics behind Figure [l] and to serve as a comparison for our subsequent 
analysis we have chosen four representative points from the parameter space for deeper investigation. 
Both points A and B fall within the cosmologically preferred region. Point B is in the coannihilation 
'tail,' so we would expect most of the annihilation cross section to come from coannihilation. As we 
can see in Figure [2] the annihilation cross section of point B is indeed dominated by coannihilation. 
We can also see the importance of t-channel sfermion exchange to x°X° ~> ff m points A and 
B, which is due to the universal scalar mass being relatively light so that this channel is open. 
Referring back to Figure [l], both points C and D are in regions where there is too much relic 
density. Figure || shows that, indeed, the annihilation is too inefficient to eliminate enough dark 
matter. The annihilation channels to two fermions still dominate but now sfermion exchange is too 
suppressed to provide the critical annihilation rate (indicated by the dashed line). As with points 
C and D, most of the parameter space of the cMSSM is experimentally excluded since there is no 
efficient way of depleting the relic density. 

It is a generic result of the low tan (3 cMSSM scenario that, excluding stau coannihilation, 
cosmological viability depends almost solely on the t-channel exchange of sfermions. This strict 
dependence causes most of the presumed parameter space to be experimentally excluded, leaving 
the dark matter prospects of the cMSSM in serious jeopardy. 

2 Nonuniversal Gaugino Masses 

The reason for the failure of the low tan f3 cMSSM to solve the dark matter problem is the low 
annihilation cross section for the neutralino. The only channel capable of providing a suitable cross 
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section is the aforementioned t-channel sfermion exchange which is only efficient in a small region 
of parameter space. If we relax the GUT relationship between the gaugino masses but still remain 
in the large \fj>\ limit (low tan/3) then we will continue to have a predominantly gaugino-like LSP 
(JVii + — 1) with the relative values of Nn and N\2 governed by the relative values of Mi and 
M 2 . Decreasing M 2 relative to M\ at the electroweak scale increases the wino content of the LSP 
until ultimately M\ > M 2 and N n ~ 0, Nu ~ 1. 

The bino component of the neutralino couples with a U(l) gauge strength whereas the wino 
component couples with the larger SU(2) gauge strength, thus enhancing its annihilation cross 
section and thereby lowering its relic density. As N12 is increased more parameter space should 
open up for dark matter until annihilation becomes too efficient in the pure wino-like limit and we 
are left with no neutralino dark matter at all Q . This is evident in Figure [| where we plot contours 
of £\h 2 as a function of scalar mass and the ratio (M 2 /Mi) at the boundary condition scale Agut- 
Allowing the gaugino masses to vary independently introduces two new degrees of freedom. We 
have chosen to vary the ratio {M 2 /M\) while fixing the value of M]/ 2 = min(-Mi,M 2 ) and M3 
at the high scale. In practice we use our choice of M1/2 to determine the value of the smaller of 
the pair (Mi,.M 2 ) at the high scale and then use the ratio (M2/M1) to determine the larger of 
(Mi,M 2 ). In Figure | we have set M 1/2 = 200 GeV and M 3 = M 1/2 . 

Particular values of (M 2 /Mi) which deviate from the universal cMSSM case allow for cosmo- 
logically interesting relic densities which are almost independent of the scalar mass. To see how 
dark matter physics changes as one departs from the cMSSM we have analyzed the six labelled 
points from Figure || in Figure |4j. Starting with the two cMSSM points, E and F, the importance 
of the process X1X1 ~~ * ff 1S again demonstrated. Point F has scalars that are too heavy for the 
two fermion final state to sufficently deplete the LSP relic density while point E has much lighter 
scalars, allowing efficient annihilation to two fermions and resulting in an appropriate amount of 
dark matter. Given the value of M]/ 2 = 200 GeV for this plot point E would lie a little to the left 
of point A in Figure [|. 

Points C and D sit at much lower values of (M 2 /Mi), resulting in a lightest chargino that is 
much more degenerate with the lightest neutralino. This enhances the importance of coannihi- 
lation channels, leading them to dominate the annihilation cross section. The main channels for 
coannihilation are x^Xi ~* //'j making up the left-hand coannihilation column in Figure |], and 
X^Xi ~~ W^Z, which is the main contributor to the right-hand coannihilation column. Chargino- 
neutralino coannihilation has become so efficient here that the relic density is now not enough to 
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account for observations. It should be noted that point D gives almost the cosmologically preferred 
relic density - its relic density is higher than that of point C for two reasons. First, the mass of 
the lightest neutralino drops slightly in going from point D (~ 122 GeV) to point C (~ 112 GeV). 
Second, the lower scalar mass at point C allows t-channel exchange of scalar fermions to go unsup- 
pressed. This is important in one of the diagrams contributing to X Xi ~^ //'• The wino content 
of the lightest neutralino is also increased by lowering (M2 jM\ ) , causing the standard annihilation 
channel to two fermions to become relatively unimportant. This increase in efficiency continues 
through points A and B, now making the neutralino relic density cosmologically irrelevant. 

The irrelevance of scalar masses above 1 TeV can be simply understood. Above 1 TeV the 
t-channel scalar exchange contribution to X1X1 ~~ > ff is suppressed due to the scalar mass. With 
enhanced wino content, however, channels that were previously suppressed are now more efficient, 
such as X1X1 ~~ * W + W~ . More importantly, we can also see from Figure ^ that much of the 
parameter space in the {Mo, (M2/M1)} plane is not ruled out: an Q x h? < 0.1 is not experimentally 
excluded, it just does not completely explain all of the needed dark matter. 



Table 1: Superpartner and Higgs mass constraints imposed [23, 24] 



Gluino Mass 


Ulg 


> 


190 


GeV 


Lightest Neutralino Mass 




> 


32.5 


GeV 


Lightest Chargino Mass 




> 


75 


GeV 


Lightest Squark Masses 


rrig 


> 


90 


GeV 


Lightest Slepton Masses 


mj 


> 


87 


GeV 


Light Higgs Mass 


m h 


> 


95.3 


GeV 


Pseudoscalar Higgs Mass 


m A 


> 


84.1 


GeV 


Charged Higgs Mass 


m H ± 


> 


69.0 


GeV 



Figure ||| examines the effect of changing the boundary scale value of M3 on the cosmologically 
preferred parameter space of Figure ||. In Figure || and subsequent plots we impose the constraint 
that the LSP be electrically neutral and that the resulting spectrum at the electroweak scale satisfies 
the search limits of Table |l]. The shaded region in panel (A) of Figure [5| is excluded by the gluino 
mass bound while the shaded region in panels (C) and (D) are excluded by the constraint on the 
stau mass. As the value of the gluino mass M3 is increased relative to the other gaugino masses 
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the cosmologically preferred range of the ratio (M2/M1) moves to slightly higher values. 

It is apparent from Figure ^ that the crucial variable in the determination of the LSP relic density 
is the value of the ratio (M^/Mi) at the electroweak scale. The region of preferred relic density 
0.1 < f^ x h 2 < 0.3 consistently asymptotes to the region between the values (M 2 /-Mi) low = 1.15 
and (M2/Mi) low = 1.25 independent of the value of M3, with the influence of the universal scalar 
mass Mq most pronounced for small values of M3. Most of the reason for this behavior is the 
composition of the low-scale squark masses. For large values of M3 the RG evolution of the squark 
masses is dominated by M3 which drives scalar masses to higher values and further suppresses 
the t-channel slepton and squark exchange diagrams. This causes the asymptotic approach to 
scalar-mass-independence to saturate for much lower values of Mq than in the low M3 case, where 
the value of squark and slepton masses is largely independent of M3 and merely a function of the 
boundary value Mq at the high scale. 

These results are robust under changes in the relative sign between the soft supersymmetry- 
breaking values of Mi )2 and M3, as well as changes in the overall gaugino mass scale Mi/ 2 , as is 
demonstrated in Figure 0. ^ For low values of M3 there is little difference between the positive and 
negative values, though for high values of M3 the entire plot moves from right to left when the 
sign is reversed. Nonetheless, the cosmologically preferred region falls between the same values of 
(M2/M\) low (denoted by dashed lines), regardless of the sign of M3: it is only the preferred region 
of (M2/Mi) high that changes with the sign flip. Both the effects of the magnitude of M3 as well as 
its relative sign can be understood from the effect M3 has on the running of M2 and M\, starting 
at two loops. The two loop running of the gaugino masses, in the conventions of [13], is partially 
given by 

where B^ b ' is a matrix of positive entries. Therefore the higher the value of \M^\ the greater the 
impact on the gaugino masses M\ and M%. Furthermore, this effect is felt more strongly by the 
SU(2) gaugino mass than the U(l) gaugino mass. Thus for a given value of (M 2 /Mi) high changing 
the sign of M3 drives the value of M2 higher at the electroweak scale to a greater degree than it 
does Mi, resulting in a higher value of (M2/Mi) low . This in turn leads to an increased relic density 
as can be seen by comparing the right and left sets of panels in Figure |7[ 

4 For the effect of changing the relative sign between Mi and M2 see H. 
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We have seen that relaxing the GUT constraint on the gaugino masses allows for significant 
improvement in the dark matter arena. This relaxation only requires a slight increase in the wino 
content of the LSP on the order of 0.1 to 5% (see Figure ||): the LSP is still predominantly bino-like 
and is not in the unappealling wino-dominated scenario which must rely on other mechanisms to 
generate supersymmetric dark matter pq| . The observations made in this section indicate that 
models which allow control over (M2/M1) at the boundary scale may be more suitable to providing 
supersymmetric dark matter than the unified cMSSM paradigm. In fact, requiring a cosmologically 
relevant relic LSP density may in turn shed light on the nature of physics at the GUT scale in models 
of supersymmetry breaking. We will carry out an example of just such an investigation in the next 
section on a class of supergravity models derived from heterotic string theory. 



3 BGW Model 

In this section we turn our attention to a class of heterotic string-derived supergravity theories 
that invoke gaugino condensation in a hidden sector to break supersymmetry. The framework for 
this model was first put forward by Binetruy, Gaillard and Wu J?], ||| and its phenomenology was 



considered in subsequent papers [|9j, |10J. As a supergravity model with a unification scale, many 
of the typical results of mSUGRA continue to hold - in particular the few number of parameters 
necessary to determine the low-energy spectrum. However, a newly-emphasized contribution to the 
gaugino masses resulting from the superconformal anomaly gives a correction to the standard 



gaugino mass unification that has been investigated recently [g7], 28]. Thus in this model one is 
able to determine the ratio (M2/M1) as a function of the parameters of the hidden sector. 

The soft supersymmetry-breaking gaugino masses M a in the BGW model are determined at the 
scale of gaugino condensation (typically of order A con( i ~ 10 14 GeV). They are proportional to the 
gravitino mass M 3 / 2 and depend on the value of the beta-function coefficient of the condensing gauge 
group(s) of the hidden sector. In practice the soft supersymmetry-breaking terms are dominated 
by the condensing group with the largest beta- function coefficient, which we label 

M M = ^^(^±M- 3 K)M 3/2 , (7) 

with 

1 / 1 \ 

(8) 
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Here C+ and are the quadratic Casimirs of the adjoint and matter representations $ , respec- 
tively, under the largest condensing gauge group. In equation (^) i is the dilaton field whose vacuum 
expectation value determines the string coupling constant at the string scale: < i >= g^ x /2. Here 
we have assumed the standard string coupling suggested by GUT scale unification: g 2 str = 0.5. 
Finally, the constants b a and b' a in (Q) depend on the group theory parameters of the observable 
sector gauge groups. For example the SU(2) gaugino mass depends on the constants: 

Csu(2) ~ g C SU(2)J 
b 2 = ^2 \CSU(2) -J2CSU(2))> ( 9 ) 

with similar definitions for the SU(3) and U(l) gaugino masses. Thus the key variable (M2/M1) 
depends on the value of b + : 

M 2 (A cond ) _ g\ (A cond ) (1 + b' 2 t) - (b 2 /b + ) (1 + b + l) 

M!(A cond ) g\ (A^) (1 + b[£) - (h/b + ) (1 + b + £) { ) 

Figure || shows contours of constant LSP relic density for the BGW model in the {6 + ,Mo} 
plane, where Mq is the usual universal scalar mass whose value is given by Mq = M 3 / 2 in the model 
we will consider here. While the axes of Figure |8| are very similar to those of Figures ||, |5| and || 
there are some notable differences in the BGW model. The gluino mass parameter M3 relates to 
M2 and M\ through an identical relationship to Equation (|h]) and therefore changes with 6+. In 
the previous figures M3 was held constant at the high scale within a single plot, but in Figure § 
the ratio M3/M1 at the condensation scale varies from 0.2 at 6+ = 0.02 to 0.8 at b + = 0.09. 
Nevertheless, there is still a region of viable dark matter largely independent of the universal scalar 
mass, as in the general nonuniversal cases studied in Section ||, for the same reasons: a smaller 
value of (M2/M1) for lower 6+ results in higher wino content as well as more degeneracy between 
the lightest neutralino and chargino, resulting in conannihilation. 

In the left plot of Figure [8| it is evident that this is not a result of the masses being tuned to 
sit on a pole. The Higgs pole, given by the locus of points for which 2m x o = m/j, is indicated by 
the uppermost dotted line. More important is the VF-pole, denoted by the second lower dotted 
line in the left plot, where a neutralino and chargino go to an on-shell VF-boson, severely warping 
the lower part of the plot. However, both of these resonant regions are excluded experimentally 
by the criteria of Table [l] as indicated in the right plot of Figure |8| by the shaded region. The key 
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constraints include the gluino mass (given by the dashed line) and the chargino mass (given by two 
parallel solid curves representing a chargino mass of 75 and 90 GeV from bottom to top). 

As with the previous cases the additonal physics contained in this plot is easily seen through a 
few representative points given in more detail in Figure H Starting from the top right in Figure |8|, 
point F sits at far too high a value of Mq for the normal cMSSM annihilation channels to be 
effective. Its high value of b + also gives a high value of (M2/Mi) low (between 1.3 and 1.4), so the 
wino content of the LSP is low. The dominant channels are neutralino-chargino coannihilation but 
this is not sufficient to deplete the relic density to acceptable levels. Point E lies exactly on the 
pole for two neutralinos going to an on-shell Higgs which then naturally decays to two fermions, 
making this the dominant final state. A lowered sfermion mass scale also allows coannihilation to 
two fermions to increase. Nevertheless, the net effect is still too small to bring the relic density 
down far enough. 

Point D is in the region where one would not expect much annihilation to fermions, but for this 
b + value the neutralino and chargino are becoming more degenerate, increasing coannihilation and 
bringing the relic density down towards the cosmologically preferred region. Additionally, points 
D and F are the only two points which kinematically allow x^X° ~~ * W^Z. Once this channel is 
open it is the main determining factor in the relic density. 

Points B and C both lie near the region where the masses of the lightest chargino and the LSP 
add up to exactly the mass of the charged VF-boson. This enhances the efficiency of most channels 
of chargino-neutralino coannihilation, resulting in a relic density that is now a little too low to 
account for astrophysical observations. For point A, by contrast, the particles are off-shell so these 
processes are too inefficient and the relic density is too high. Note that for points A, B and C the 
value of the lightest chargino mass is below the experimental limit so these points are excluded. 



Figure 1C shows how the parameters of the BGW model determine wino content and the ratio 
(M2/M\)-y . As in the case of Figures [5| and ^ from Section ^, cosmological observations favor 
a mild wino content of 0.1 to 5% and single out the region 1.15 < {Mi/M\)^ < 1.25. The 
correspondence between the value of b + and (M2/Mi) high is clear from the comparison of the right 
panel in Figure 1C and those of Figures |5| and ^, in particular panels (A) for lower values of 6+ and 
(B) for higher values. 

To see the discriminatory power that cosmological considerations can have on model building 
we now look more deeply into the role of the hidden sector configuration in determining the pattern 
of soft supersymmetry-breaking terms in this class of supergravity models. The (nondynamical) 
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gaugino condensates in the hidden sector are represented by dimension three chiral superfields 
U a — Tr(W Q W Q ) a where a labels the condensing groups of the hidden sector: <5hid = Yl a Qa- 

The superpotential for these low-energy effective degrees of freedom is that of Veneziano and 
Yankielovich [E9(l 



+ h.c, 



(11) 



which is here written in the chiral U(l) superspace formalism The lagrangian involves the 

gauge condensates U a , the complete Kahler potential K, and any possible gauge-invariant matter 



where a runs from 1 to the number 



condensates described by chiral superfields IP ~ ]J A (3> J 
of condensates N c . The coeffecients b' a , b% are determined by demanding the correct transformation 
properties of the expression in ( |i~l| ) under chiral and conformal transformations (?], 3C] and yield 
the following relations: 

ba - b' a + TK = ^ 2 U ~ lT C a) ■ Y»i P a = 

a °" V A / a,A A ' 

which are equivalent to those of @ and (||). The matter condensate superpotential is taken to 
be V7[(n p< *),T] = J2 a c a W a (T) (IP) Pa , where T represents one of the three untwisted moduli 
chiral superfields T 1 which parameterize the size of the compactified space. The coefficient c a is a 
Yukawa coefficient from the underlying theory which we presume to be O (1). Finally, we require 
on dimensional grounds that p a J2a n a = 3 for all values of a. 

Armed with these elements of the lagrangian the equations of motion for the nondynamical 
condensate superfields can be solved f7j 3 and the condensation scale and gravitino mass deter- 
mined. If we define the lowest components of the chiral superfields as u a = U a \ 9= Q =0 = p a £ lu)a and 
tj = Tj\ 9= q =0 then the condensate value is given by 

Pa 2 = e- 2 ie K e~^Z e ££, ln L A/\ ^ jj |^ /4Ca |- 2 f , (13) 



where ^(t 1 ) is the Dedekind function. To disentangle the complexity of (13) it is convenient to 
assume that all of the matter in the hidden sector which transforms under a given subgroup Q a is of 
the same representation, such as the fundamental representation, and then make the simultaneous 
variable redefinition 



cff 



N, 



(14) 
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In the above equation 6J, ep is proportional to the quadratic Casimir operator for the matter fields 
in the common representation. 

From a determination of the condensate value p the supersymmetry-breaking scale can be 
found by solving for the gravitino mass, given by M 3 / 2 = \ (\J2a^aUa\) 01, though in practice we 
will replace the summation with the condensing group with the largest beta- function coefficient: 
M 3/2 = ±6 + (|u+|). Now for given values of (c a ) eS the gravitino mass can be plotted in the 
(fr+) cff j plane, as in Figure 11, where curve (a) is a contour of M 3 / 2 = 100 GeV for (c a ) cS = 10 
and curve (b) is a contour of M 3/ / 2 = 10 TeV for (c a ) cS = 0.1. The shaded region between 
these curves can then be thought of as the phenomenologically preferred region of hidden sector 
configuration space. 

Upon Zn orbifold compactification of the heterotic string the £"8 gauge group of the hidden 
sector is presumed to break to some subgroup (s) of Eg. For each such subgroup the equations 
in ( |l2| ) define a line in the {ba) e s} plane which we have displayed in Figure 11. We then 
sampled 25,000 combinations of (6") off , (c Q ) eff } which give rise to gravitino masses between 
100 GeV and 10 TeV and which yield a particle spectrum consistent with the bounds in Table |]. In 
Figure [Tl] we display those combinations which implied a relic density in the range 0.1 < Q x h 2 < 0.3 
(fine points), as well as the slightly higher range 0.3 < f2 x h 2 < 1.0 (coarse points). 

Figure |ll| clearly favors a very specific region of hidden sector parameter space with a preferred 
value of 6+ in the neighborhood of 6 + = 0.036 and a corresponding range in (&" ) eff of 0.2 < (b") cS < 
0.6, which points towards a large condensing group such as 50(12), 50(14), 50(16), E§ or Ey. A 
typical matter configuration for the hidden sector would be represented in Figure [II] by a point on 
one of the gauge group lines. The number of possible configurations consistent with a given choice 
of {a s tr, (c a ) eS } and supersymmetry-breaking scale M 3 / 2 is quite restricted. For example, if we 
ask for a hidden sector configuration charged under the Eq gauge group for which Ce 6 = 12 and 
0|^ nd = 3, and require that our matter condensates be gauge invariant so that fundamentals must 
come in groups of three, then from (|l^) the only combination that falls in the preferred region 



of Figure |ll| is Vf unc j = 9. This combination is notable in that it was shown in ||10[ to possess 
many desirable phenomenological features. A similar analysis for the other allowed gauge groups 
leaves only a handful of possible hidden sector configurations, summarized in Table |2j, where we 
have included some examples with various hidden sector effective Yukawa couplings (c a ) eS and the 
implied values of M 3/ / 2 and ^h 2 - ^ s ^ s evident from the table and from Figure [Tl], using the 
dark matter constraint on LSP relic densities is a very powerful tool in restricting the high energy 
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physics of the underlying theory. 



Table 2: Gauge group Casimirs and allowed condensate numbers. 



Gauge group 


C a 


£<fund 


b a 






( c a)eff 


M 3/2 (GeV) 


n x h 2 


E 6 


12 


3 


0.038 


0.23 


9 


3.8 


5967 


0.633 


50(16) 


14 


1 


0.034 


0.29 


34 


2.7 


7011 


0.194 


SO (14) 


12 


1 


0.034 


0.24 


28 


4.4 


3383 


0.069 


SO (12) 


10 


1 


0.034 


0.19 


22 


6.3 


1438 


0.076 



Conclusion 

The prospects for cMSSM dark matter are rapidly diminshing, barring a curious conspiracy between 
Mq and Mi/2- This is due to the inefficient annihilation of a dominantly bino-like LSP. Departure 
from the standard cMSSM GUT relation allows values of (M^/Mx) that accomodate small admix- 
tures of wino content for the LSP. Lowering this ratio at the electroweak scale increases the LSP 
annihilation efficiency by virtue of its higher wino content and the tightening degeneracy between 
the lightest chargino and the LSP, resulting in increased coannihilation. Ranges of (M2/Mi) low 
exist with 0.1 < Qh 2 < 0.3 and where the value of Mq is restricted to be anything above 1 TeV - 
quite in contrast to the very light scalars required in the standard cMSSM case. 

The requirement of cosmologically interesting relic densities, or at least the demand that ^ x h 2 < 
1, can be a powerful constraint on models with nonuniversal gaugino masses which is often quite 
complementary to the constraints arising from direct search limits for superpartners. As an example 
we investigated the BGW model of gaugino condensation derived from heterotic string theory where 
the number of possible hidden sector gauge groups and matter configurations could be restricted 
to a very small number. Similar analyses on models with small deviations from universality should 
prove equally fruitful. While relic densities of supersymmetric particles that were once in thermal 
equilibrium need not be the explanation for the missing nonbaryonic mass in the universe,^ it is 
nevertheless one of the most compelling aspects of low-energy supersymmetric phenomenology and 

5 In for example, nonthermal mechanisms are used to provide adequate relic densities in the case of the highly 
wino-like LSP characteristic of the standard anomaly-mediated supersymmetry breaking scenario. 
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promises to remain so even in scenarios with heavy squarks and sleptons. 
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Figure 1: Preferred Dark Matter Region for cMSSM. Contours of Q x h 2 of 0.1 (bottom-most contour), 0.3 
and 1.0 (top-most contour) are given. The shaded region is ruled out by virtue of having the stau as the LSP. The 
Higgs pole region and stau coannihilation tail are clearly discernible. We have also added contours of constant Higgs 
mass for rrih = 100 GeV and mi = 115 GeV. The four labeled points are examined in Figure 0. 
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Figure 2: Annihilation Cross Sections for Selected Points from Figure These graphs detail the com- 
position of the neutralino depletion cross section. The total depletion cross section is on the far left. The next two 
columns divide the normal annihilation channels by final state into two fermions or all other annihilation channels. 
The final column is the sum of all coannihilation channels. The total relic density is given at the top of each plot and 
the dashed horizontal line illustrates the ideal total depletion cross section for an S7 x /i 2 = 0.2. 
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Figure 3: Preferred Dark Matter Region for Nonuniversal Gaugino Masses. Contours of Q x h 2 of 0.01, 
0.1, 0.3, 1.0 and 10.0 from left to right, respectively, are given as a function of the ratio of SU(2) to U(l) gaugino 
masses M2/M1 at the high scale. The cMSSM is recovered where the two masses are equal at the high scale, as has 
been indicated by the dashed line. The six labeled points are examined in Figure H. 
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Figure 4: Annihilation Cross Sections for Selected Points from Figure g. These graphs are identical in 
layout to those of Figure ^, except now the coannihilation channels are split into two columns: two-fermion final 
states and everything else. 
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Figure 5: Preferred Dark Matter Region with Nonuniversal Gaugino Masses for Various Ms Values 

I. Contours of constant relic density are given by the solid lines for f2 x h 2 = 0.01, 0.1, 0.3, 1, and 10 from left to 
right. The dotted lines are curves of constant wino content for Nl 2 = 0.25, 0.04, 0.01 and 0.025 from left to right. 
The value of the high-scale gluino mass M3 is given in terms of the universal gaugino mass M 1 / 2 at the top of each 
plot. This plot uses a value of M 1 / 2 = 200 GeV except for panel (A) where M^/ 2 = 250 GeV. The shaded region is 
excluded by the constraints of Table |XJ. 
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Figure 6: Preferred Dark Matter Region with Nonuniversal Gaugino Masses for Various A/3 Values 

II. Contours of constant relic density are given by the solid lines for f2 x h 2 = 0.01, 0.1, 0.3, and 1.0 from left to right. 
The value of the ratio {M^/Mi)^ is indicated by the dashed lines for the values (M2/M1), = 1.15, 1.25 and 1.50 
from left to right. The shaded regions are ruled out because of a stau LSP. 
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Figure 7: Preferred Dark Matter Region with Nonuniversal Gaugino Masses for Differing Signs of 

M3. The solid lines are lines of constant f2 x h 2 and the dashed lines represent values of (M2/M\) low = 1.15, 1.25 and 
1.50 as in Figure []. The shaded regions are excluded by virtue of having a stau for the LSP. We have indicated the 
Higgs mass constraint in panels (A) and (B). Here M x / 2 is taken to be 500 GeV in all panels. 
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Figure 8: Preferred Region in the BGW Model. In the left panel contours of constant relic density are given 
as a function of Mo and fe+ by the solid lines while the Higgs resonance and W resonance are indicated by the upper 
and lower dotted lines, respectively. Moving outward from the lower dotted line are contours of f2 x h 2 = 0.01, 0.1, 
0.3, 1.0 and 10. In the right panel these contours are reproduced with experimental constraints from Table |lj The 
shaded region is excluded: the dashed curve represents a 190 GeV gluino mass while the two parallel solid curves 
represent a 75 GeV and a 90 GeV chargino mass from bottom to top. The labeled points are examined in detail in 
Figure ^[ 
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Figure 9: Annihilation Cross Sections for Selected Points From Figure g. These graphs are identical in 
nature to those of Figure 0. 



26 



M o (GeV) 
6000 



5000 



4000 



3000 



2000 



1000 




M o (GeV) 
6000 



5000 



4000 



3000 



2000 



1000 




0.02 0.03 0.04 0.05 0.06 0.07 0.08 

b + 



0.09 



0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 



b + 



Figure 10: Preferred Region in the BGW Model. Contours of constant relic density are given as in Figure |^ 
by the solid lines for f2 x h 2 = 0.1, 0.3, and 1.0 only. The left plot gives dotted lines of constant wino content (25%, 
4%, 1%, 0.25% from left to right). The right plot gives dashed lines of constant ratio (M 2 /Mi) low (1.15, 1.25, 1.5 
from left to right). 
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Figure 11: Preferred Dark Matter Region in Hidden Sector Configuration Space. This plot illustrates 
the dark matter parameter space in terms of the gauge group and matter content parameters of the hidden sector. 
The fine points on the left have the preferred value 0.1 < fl x h 2 < 0.3 and the coarse points have 0.3 < f2 x h 2 < 1.0. 
The swath bounded by lines (a) and (b) is the region in which the 0.1 < (c Q ) eff < 10 and the gravitino mass is 
between 100 GeV and 10 TeV. The dotted lines are the possible combination of gauge parameters for different hidden 
sector gauge groups. 
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